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Аналітичне продовження сили тяжіння в регіональному масштабі визначає розв'язок 
нелінійної межової задачі Длексідзе для рівняння сили тяжіння. Оцінено похибку її заміни 
на відповідну задачу для рівняння Лапласа. Точність розв'язку істотно залежить від міри 
обумовленості задачі, наближень напрямку шуканого градієнта, параметрів чисельного ме- 
тоду, напрямку зовнішньої нормалі, вздовж якої обчислюють похідні у граничних умовах. 


Ехасі дгаміїу ітап5Ї0гтп5 їп пе гедіопа! 5їидїез аге деїегтіпеа Бу їпе 50ГПийоп ої попіїпеаг 
Боипдагу-уаїше Діехідге ргобіета Їог пе дугауіїу роїепіа!. Тре егтогя ої ії5 зар5ійийоп Бу Ше 
арргоргіаїе ргобіеті Їог бе ІТ аріасе'8 едцайоп Бауе Бееп е5бітаїеа. ТРре асситасу ої ії8 50Їа- 
оп ехігетеї!у дерепдз оп Ше ргобіет'я сопаїйопіпд, арргохігайопя ої бе гедцігед дтаді- 
епі дігесіїоп, рагатеїег5 ої питетіса! саїсиїйц5 апа дїігесіїоп ої їБе ехіегпа погта!, а|опд 
місь бе ргорег дегіуайуве5 аге саїсціаїеай. 


Трансформации силь тяжести. Задача восстановления потенциала силь тяжести в боль- 
ших областях возникла еще в трудах Клеро о фигурах равновесия в виде сфероидов, которьтве 
обосновали гравиметрию как науку. Последующие работьт Стокса, Пуанкаре и Молоденского 
указали путь ее решения в виде смешанной граничной задачи для уравнения Лапласа. Из-за 
приближенности граничньх условий зта задача гарантирует удовлетворительную точность вь- 
числений потенциала только в локальньтх областях "малой мерьт" |Якимчик, 2001). 

Другой путь вьшчислений потенциала основан на том, что, согласно работе |Юнг, 1936), 
аномалий сильг тяжести считают в первом приближеними значениями " вертикальной " произ- 
водной зтой аномалий, т. е. реальнье значения сильт тяжести замененьг гармонической функ- 
цией У,. Ото позволило на базисе теории гармонических функций возвести всю теоретичес- 
кую надстройку современной гравиметриий. Но погрешность такой заменьі, как указано в ра- 
боте |Якимчик, 2001| и ниже, растет с размерами области сосредоточения аномальньтх масс. 
Вследствие сказанного для надежного восстановления потенциала при региональньгх иссле- 
дованиях следуєт воспользоваться иньми методами, например |Дубовенко, 2009). 

Гравиметрическая сьемка определяет набор значений модуля градиента потенциала и си- 
ль тяжести (МГПСТ) 2| фа ятад И | 57М (5) на некоторой замкнутой дневной поверхности 5, 
которая целиком охватьваєт тяготеющие массьгт в области С поиска трансформант. Прямая 
его трансформация не имеет отношения к линейньгм граничньгм задачам математической фи- 
зики и является нетривиальной задачей из-за неизвестного направления вектора 2 5: При 
приближенном решений зтой задачи ее граничнь6е даннье можно задавать на незамкнутой, 
но довольно большой поверхности |Алексидзе, 1965. 

Трансформацию значений гравитационного поля определяют путем решения граничной 
задачи непосредственно для нелинейного дифференциального оператора А (є) силь тяжести 
в области задания є|х ІЧернькй, 1970): 


А(д) «44 (хі ) - а" (9, хі )9(х )зу(х), хеб, і 21,3, (1) 
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3 З 3 2 
9-8 5 
где а? ( 9,Хі )з р со5 Ці Д со5 0 з У ЬУ зх Й -- геометрический фактор поля 
ізі різі їзі хі 9 
тяжести, У(Х| ) - - АЛУ | ягад є (Хі )| -- заданная функция распределения масс, й; -- Угол 


между направлением вектора 2 5 Й і-й координатной осью. Зто уравнение учитьваєт на- 
правлениєе сильг тяжести и сложно для вьічислений. Если вектор 2 5 В любой точке прост- 
ранства направлен в одну и ту же точку Р ( хо ),і 21,3 (вдоль нормали к поверхности зл- 
липсоида), уравнение (1) упрощаєтся |Алексидзе, Гелашвили, Картвелишвили, 19721: 





лу (хі )-29(х)/ г? зику), ге ДУ (жо). о) 


Практические способь перерасчета значений сильт тяжести -- постановка и решение гра- 
ничной задачи непосредственно для значений МГПСТ |Алексидзе, 1985| либо решение нели- 
нейной внешней граничной задачи Алексидзе для потенциала сильгт тяжести: 





Линеаризация задачи (3) возможна в предположений касательно направления МГПСТ на 
поверхности 5, которое вносит значительнье погрешности, а определение направления силь 
тяжести -- громоздкая задача, как и определение МГПСТ. Для внешней тяготеющей точки 
получено уравнение |Алексидзе и др., 1972| 


3 


9Д со8 53 9 
ізі 


д со8 їн 
що ЗАдсоваї-0, Ііз-і1,3, 
дХі 


определяющее модуль сильг тяжести є по известньїм распределениям углов (і в простран- 
стве С (или косинусь углов а; по известньжм значениям 2 в С с неопределенностью со58 а; на 
бесконечности). 

Вопрос разрешимости и единственности задачи (3) для потенциала сильг тяжести принци- 
пиален в гравиметрии: отвечаєт, возможно ли вьісокоточное однозначное определениєе силь 
тяжести во внешнем пространстве по точньтм абсолютнькм значениям сильг тяжести на замк- 
нутой поверхности 5, целиком охватьтвающей все массь(к, и какие нужнь минимальньтее допол- 
нительнье даннье. Вьічисление сильт тяжести во внешнем пространстве при условиях одно - 
временного определения на поверхности 5 ее абсолютньгх значений и направления неслож - 
но, но дорого и делаєтся редко. 

Решают задачу (3) путем последовательньжх приближений. В плоском случаєе (граничньге 
значения м (5) й функция и зависят от двух переменньгх) задача определения гармоничной 
функции, модуль градиента которой задан на замкнутом контуре, 


лц(х|)-0, хієб, иц| 20, 
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не имеет единственного решения |Алексидзе, 19851. Если 8 (2) и (Хі ) я ЇМ (Х; ), і з 12 
произвольная голоморфная в С функция, действительная часть Ве /(2) которой совпадаєт с ре- 
шением граничной задачи (4), ее мнимая часть Пі / (2) также является решением |Алексидзе, 
1985 | задачи (4). Отношение Коши -- Римана К| : К» - - 1 - - условие перпендикулярности 








ди ди ду ду 
градиентов и, у, где К| з і8 4 З Зх СР Ка з ів аз ах У угловьте козффи- 
2 1 2 1 


циентьг касательньгх вдоль градиентов функций и и у (Є|, (2 - - угльт между положительньтм 
направлением оси и и и у). 


ди 
В конечной и-окрестности направления / градиента єгад - --- одного из решений гра- 


ді 
ничной задачи (4) не существуєт другого решения задачи (4), т. е. нет двух решений задачи 
(4), направления градиентов которвбіх составляли бьі угол меньше й (угол между градиентами 
двух решений (4) - 1 / 2, решения, различающиеся постоянной, -  тождественньт). Проверка 
зтого предположения осуществлена путем численного решения методом | Алексидзе, 1978| 
граничной задачи (4) с помощью итераций 





-1) ди ПО 
дик (хі)є0, У с | ; це мів), і еАї, 2, 5) 
ї а 
где козффициентьт 
2 
би ан 
Ко 
са (и )- дхі 2. дхі 


направляющих косинусов полученьт нормировкой граничньгкх условий задачи (4), й Ю оз бре 
приближение ее решения. 

Погрешности трансформаций. Погрешность є заменьт уравнения (2) на уравнение Лап- 
ласа 9 (Хі ) - 0 -- разность ренанни на сфере С радиуса а граничньх задач для уравнений 
силь тяжести Д9 (Х|;)- «а СХ ур :0,хєе С, и Лапласа 4 (Х; )-0,хє С, с одинаковьім 
граничньтм условиєм 5 | 5 З у (5). вактов є | 5 направлен в Одну и ту же точку О (5, п, 0). 
Погрешность удовлетворяєет граничной задаче Дє - 24 (Хі )/ г? хе с, с однородньми гранич- 
ньми данньми є Р - 0, причем и 3 є Х у, Если начало координат совпадаєет с центром сферкбі 
С, то в точке Х; з 0, і 21,3, имеем 





а? / 6чє ка?ч, / 6, (6) 


и при условий, что область С касается поверхности Земли, нормальное гравитационное поле 
будет равно 





1 2 
/ 2.981 
у, 2 тах 20(х|)/ г2 «2 тах 9(х; )/ піп Г со ЗВАНІ 
іє б іє б / |хіеб (6360)? 
заз іду бабо аіь БЬЛН задо за 2 -(978000-300а). 
хіє б хіє б / хіє б (6380-- а)? 


Отсюда, согласно (б), при а- 100 км получаем 76 «є 481 мГал. В перерасчете сильт тяжести в 
свободном пространстве по формуле Пуассона погрешность может достичь 100 мГал на боль- 
ших площадях земной поверхности. 
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2 / 2 


В плоском случае У2 а?- У хо /А«еху а?- У хо 1/4 при хі2 0,213, 


і ші / ізі1 / 
даєт 514є 454 мГал. | | 

Согласно (б), є 9-0 4 0, позтому действительное значение сильт тяжести всегда меньше, 
чем полученное в предположений ее гармоничности (на зто указьваєт и де20 -- условие во- 
гнутости поверхности). 

Если на сфере С величина 2/г 7 постоянная, то погрешность є; является решением гра- 


2. 4/кЗ 


0 Ї ге б, где Го -- среднее значение г в сфере С. 


ничной задачи Де| - 29 п/ г 


3 
1 2 
Вьгражениє ГУ дно ьЗ а но 0,03 мажорируєт є,, позтому є; «2, мГал, а приах 
ба 6370 
із1 


«100 км -- Є; 4 1 мГал. Таким образом, членьт при " в уравненим (2) почти постояннь. Срав- 
нение соответствующих погрешностей указьтгваєт, что уравнение сильт тяжести лучше аппрок- 
симировать уравнением (2), чем уравнением Лапласа. 

Уравнениєе сильт тяжести (2) при / - Го бьштрее сходится, а при малом а /Гд погрешность 


3 / 
зтой заменьтг є2 - 20 Х3 а"- ру Хі / 5 о имеет зкстремумь Х; 2 Хо з 0, Хз с- Жа/ / 3. 
інаі | 


На границе 5 аномальной области С є» є 0, т.е. тах є» | -4да ці 15 3 п откуда приагз- 
-100 км є» «0,5 мГал. В зтих оценках погрешности заменьг уравнения (2) уравнением Лапла- 
са предположено, что на замкнутой границе 5 известньт точнье значения силь тяжести; при- 
ближеннье значения 2 |, несколько увеличат величину погрешностей. 

В практике гравиметриим при перерасчете во внешнем пространстве решают не внутрен- 
нюю граничную задачу, а внешнюю (ее сводят к решению внутренних граничньх задач |Але- 
ксидзе, 19651). Погрешность перерасчета во внешнем пространстве С (заменьт внешней гра- 
ничной задачи Лд (Х|; ) - 24 (Х; 5/1 - 0 на внешнюю задачу Л9| (Х; )- 0 при одинаковьіх 


С 
граничнькх условиях 2 |, , 2 У(5), 5|,, 7 0)равна є з - (з Пре іт) и достигаєт макси- 
го оо а 


мума в области С при условий є | грошо 79 (а) /4, где 2 (а) - - значение сильт тяжести на по- 
верхности сфероида. Предположениєе гармоничности сильт тяжести в перерасчете на вьісоту 
6000 км дает погрешность - 245 Гал. 


д М М ( г - (со50 ) 
Радиальная компонента У, з - в 


дг 
нка)? (г2 -огосовоч С? 








а 


аномалий У (Х;)-М / б з х з (Х3 - 0) 2) от однородной сферьг массой М с центром в 


- С 
точке Р (0, 0, б « 6375) имеет погрешность перерасчета є - - -У,, гдес -- произвольная по- 
Г 


стоянная. В точках г-а и 0-0 погрешность равна є; - М |-а / (а - уз г--У, |, аее макси- 
мум при любом фиксированном г существуєт в 0-0, возрастая от - 2е(г - 0) до т со (Х3 - С), 
спадает к отрицательному максимуму в точке г - г (0) и равен нулю в точке г - а. Если масса 
М -100(а-0) 2 аномалия в ближней точке поверхности (0, 0, а) равна 100 мГал; при а- 6371, 
2-6361 и г - 6381, є; - 8 мГал. 

Итак, трансформацию наблюденного или нормального поля силь тяжести во внешнем 
пространстве можно осуществлять путем решения методом |Алексидзе, 1978| граничной зада- 
чи для уравнения сильт тяжести (2), конкретнье операторкеі которого для разньх референц- 
поверхностей |Алексидзе, 1985| в случае земного зллипсоида вращения (с достаточно мальнюм 
сжатием) не генерируют существенно различньгх результатов близ Земли. 
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Необходимь6е в численном моделированиий явнье аналитические вьтражения фундамен- 
тальньх решений уравнения силь тяжести (2) полученьг в общем виде |Черньтгй, 1982|. Для 
вьщчислений важно представлять соответствующие интегральт в злементарньхх функциях. 


Перерасчет "аномалий" сильг тяжести 62 - є (5) - у (5) (разности нормального и наблюден- 
ного поля силь тяжести) имеет свой нюансь |Дубовенко, 2010). Впрочем, с точностью 0 «є 2 
« 0,04 мГал аномалий сильг тяжести на равнине являются значениями вертикальной составля- 
ющей по направлению модуля градиента. 


Перерасчет в нижнее полупространство известной плотности с (Р) в уравнений (2) несло- 
жен. В практике региональной гравиметриим при перерасчете в сторону масс нередко на по- 
верхности 5 известньт вертикальнье составнье д// / дХз аномалий, а под поверхностью ано- 
мальнью массьт расположеньтг на неизвестном расстояниий от нее и следуєет гармонически 
продолжить д// / дХз в сторону масс к ближайшим особьгм точкам поля. Задача в подобной 
постановке существенно некорректна (имеет бесконечное множество решений |Алексидзе, 
19851). 


Результатьт вьічиСслений. Численно решена |Алексидзе, 1978| нелинейная граничная за- 
дача (3) для пространственного случая с помощью итерационного процесса 





3 
дип 2 
дир з 0, | У ої (шо) 3.3 | зм" (5), (Р 
і е1 5 
где й» и Ц5-;7- пе (п - 1)-е приближение, й (йп-1) -- ій направляющий косинус гра- 


диента (п - 1)-го приближения й; - |. С целью проверки сходимости итерационного процесса 
(9) решеньг такие внешние граничнье задачи: 


Мі-о, ва 2 в 


2 2 2 
Хі Хо кХо 


мі-о, |втаи, - і - і | 9 


(хоп ху юку 0 (ху Р)? ю хо хо 








2 2 2 
2 (о кхіззк) 


(ху -ву жк хх в) ніч) 





ДР 3 р ДР 3 
(хіх -вхуч) (з ж; -вхучю) 
з з , (10) 


(ба -пуд хіх) (уж пуд хі ні) 





где й - 10, Хз з 5 км, 5 - - горизонтальная плоскость. Точное решение граничньх задач (8) и 
(10) таково: 
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й й й 


ГЛИН ЗР гори рн во офо зі 
ро бот) зак Ч (а п) кох 


а граничная задача (9) точного решения не имеет. Рис. 1 иллюстрируєт результатьгт решения 
за 5 итераций задачи (8) в 30 точках поверхности 5 на равномерной сетке с шагом ДХ; - 10 км 
в области - 255 Х» 225 км, - 252 Х| 5 15 км (нумерация узлов -- рост Х| при нарастаний Х»). 
Первое приближение направления грубое -- й (й)2 02 (ц1)-0, аз (ц|)- 1. В другом 
варианте решения взято й (й) - 0,8, 0» (ц)-0, аз (ц;)- 0,6, а градиент й, и правая 
часть граничного условия (8) совпали с точностью до є - 10 6 за б итераций. Модули точного 
и приближенного градиентов совпадают с точностью 10 7. 

Граничнье задачи (8)--(10) решали для разньх первьїх приближений и положения вспо- 
могательньх точек | Алексидзе, 1978|, определяющих фундаментальнь»е решения. В случає 
рис. 1 вспомогательнь6е точки взятьтг на плоскости Хз - 1. Число итераций в вариантах варьи- 
ровало незначительно (не более 10). Граничнье значения точньхх и приближенньтх градиен- 
тов совпадают с точностью є - 10 7", они сами совпадают с меньшей точностью из-за плохой 
обусловленности задач с производньми в граничньгх условиях (для них, в отличие от задачи 
Дирихле, несправедлив принцип максимума). Оту особенность, не связанную со сходимостью 
итераций, следуєет учитьвать в практическом моделирований. 

В определений погрешности решения задачи (8) задействованьт 36 фундаментальньх ре- 
шений, равномерно расположенньітх на плоскости Х3 - 1, граница 5 - - плоскость Х3з - 2. Ле- 
вая и правая части граничньх условий (8) после четвертой итерации совпали с точностью є - 
- 1076, относительная погрешность решения и -- 10 9, абсолютная -- 10 ; в производньх 
погрешности достигают 80 9. Увеличение числа функ- 
ций до 49 снизило погрешность решения и, производ- 
ньх -- до 3--5 9, ди / д Хз -- до 20 9, Наращивание 
их числа до 64 увеличило погрешность решения и 
производньх до 70 9, хотя правая и левая части гра- 














т Составляющая градиента 
675-423 


6254384 (| 
575 





30- / ничного условия (8) после восьми итераций совпали 
525 с точностью є - 10 2, Зтот зффект обьясняют |Алек- 
475 сидзе, 1985| неминимальностью систем фундаменталь- 
425 ньжх решений и малостью детерминанта матрицьіт ко- 
375 зффициентов для фундаментальньгїх решений. 
зав Точка Х; - Х, є Хз - 0, создающая поле в задаче 
(8), - -одна из вспомогательньх |Алексидзе, 1978|, оп- 
Яга ределяющих фундаментальньгюе решения. Зтим обь- 
225 


ясняем немалую точность определения производньх 
ж їв я Зб, хо ди/дх, по приближенному решению. Если она не 
І» 16 20 М 28 входит в состав вспомогательньсх точек, точность оп- 
ределения производньтх хуже, особенно для ди /д Хз3. 
Если направление градиента на всей границе 5 
образует с внутренней нормалью острьсй угол, внеш- 
няя граничная задача (3) решаєтся с удовлетворитель- 
ной точностью и погрешность обусловлена не плохой 
сходимостью итераций (7) (число их мало), а плохой 
обусловленностью граничной задачи с косой произ- 
водной. 
В общем случаєе, граничная задача с косой произ- 
водной вдоль направления І 





ди (5) 
Рис. 1. Решение внешней граничной Ди -0, ЕТ "М (11 
задачи Алексидзе (8) и фрагмент соост- ; 
ветствующей таблицьт даннькх (х 103). | при условий, что направление внешней нормали М 
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неортогонально направлению | (/ ни в одной точке границь С не совпадаєт с касательной) 


іє (М,1)»0, (12) 
х» є 5 
не имеет решения, отличного от постоянного |Бицадзе, 1966). В граничньх задачах (8)- (10) 


зто условие вьшолняєтся на любой конечной части плоскости 5, так как поле создано точеч- 
ньшм источником, а на всей бесконечной плоскости -- нет, ибо в бесконечно удаленной точке 
плоскости 5 (М, І) | і З 0. Для внутренней граничной задачи (11) и точечного источника 
Р вне области С условие (12) для достаточно гладкой границь невьштолнимо, так как всегда 
будет существовать точка на 5, в которой вектор с началом в Р ортогонален к нормали. Чи- 
сленнье зкспериментьт на сфере С обнаружили плохую обусловленность граничной задачи 
(11), хотя граничнье условия удовлетворяются хорошо, вьтчисленнь6е и точньг?е значения гра- 
диента и и его производньх существенно разньге. 

Решениєе внутренней граничной задачи (8) для единичного куба С (внешняя нормаль имеет 
шесть разньхх направлений, начало координат (0, 0, 0) расположено на расстояний 0,2 от ниж- 
ней кромки куба С, а 24 вспомогательньге точки взятьг на кубе со стороной 1,4) согласно (7) 
при о (Цц5)- а» (Цц5) - 0, дз (4) - 1 таково: после пятой итерации граничнье значения удов- 
летвореньі (є « 1 95), но приближениєе и и его производнье не имеют ничего общего с точньтм 


решением 0) «й Ю, Ж З б з Ж и его производньми; увеличение числа итераций до 24 не 


увеличило точность приближений. В табл. 1 в четьгрех точках границь (точки коллокации ) 
даньг точнье и приближенньте значения градиента и производньх (в других точках они из-за 
симметрии совпадают с одним из приведенньхх в табл. 1). Добавив 25-е фундаментальное ре- 
шение (оно совпало с точньм решениєем задачи (8)), после четвертой итерации граничньге ус- 
ловия, решение и производнье стали удовлетворять с точностью є - 10-73, Чтобьт разделить 
влияние неустойчивости граничной задачи (11) на суммарную погрешность решения гранич- 


ной задачи (8), для функции 0 - й ГИ Ж з й з Х, решена граничная задача (11). Послед- 


ние столбцьг6 табл. 1 дают решение й граничной задачи (11) в тех же точках: погрешность 
обусловлена неустойчивостью задачи (11), ане медленной сходимостью итераций (7). 
Сходимость итераций (5) в плоском случае имеет особенности: точность решения внешней 


граничной задачи (4) с граничньм условием ї (5) - 1 / М І з Хо - 1 иточньм решени- 


5 
ем йи-їп З з Х, для единичного круга с центром в начале координат при начальном 





направлений, а ш со5 (КЕ -Ф), о) з 5іп (Її 4 Ф) (7 -- центральньвй угол, Ф -- возмущение на- 
чального направления) зависит от Ф и М. В табл. 2 приведено количество итераций М до 


- (к) -5 з 
условия их остановки | Є | -|| стад - 11410 ", максимальнье погрешности решений 





в точках коллокации после каждой итерации для производньжх -- є; « 10 а между точками 


Таблица 1. Решение внутренней граничной задачи (8) и задачи с косой 
производной (11) для единичного куба 

















4 9 3 2 6 30 18 8 13 1 2 8 


142 Геофизический журнал Мо 1, Т. 33, 2011 


ОБ ОПРЕДЕЛЕНИЙ ПОГРЕШНОСТЕЙ ГРАВИМЕТРИЧЕСКИХ ТРАНСФОРМАЦИЙ 


Таблица 2. Параметрь решения плоской внешней граничной задачи (4) 
для единичного круга 





0,000003 | 0,005021 | 0,000004 | 0,850805 | 0,000005 | 0,999000 | О,О00060 | 0000004 


ди ди 02) ди 
- тах 
І 





9х 9 Хі 9 Хі 

Решение внутренней граничной задачи и при благоприятньжх условиях определяєтся с 
точностью до постоянной (позтому не сравниваем точное и приближенное решения) и знака 
(при ф- ля первая итерация удовлетворяєт граничной задаче). 

Решена задача (4) с радиусом круга г - 2 и направлением градиента ар з щ з 2/2 
(совпал с прямой Х;| 2 Х»). Зто направление образуєт с нормалью угол 907 в точках круга 


час ма В 
о 


ТВ Пре ка РА ДИ РР 
1179 УЗ но п 


коллокации в конечном решений є» е що. где 812 


, т. е. не вьшолняєется условие (12). Точное решение 


задачи и - 0,693, а граничное условие м (5) - 0,5. В приближенном решений задачи (4) (рис. 2 -- 
результат девяти итераций, где Х|, Х» - - коОрдинатьт граничньгх точек коллокации) градиен- 
ть приближенного и точного решения близки, а сами решения и их производнье существен- 
но разнье, так как приближениєе б щ о щ 12 /2 не обеспечиваєт условия (12). 

Для сходимости итераций (5) необходимо не только, чтобьт тшочное решение удовлетворяло 
условию (12), но и вьбор начальньжх приближений, для которьїх в итерациях (5) не возникали 
бь козффициенть 0; , не обеспечивающие условия (12). 

Решение по итерациям (5) внутренней граничной задачи для единичного круга с центром 





в начале координат, граничньм условием їм (5) з і М; (хр - 1) 21 (х» - 1) 2 при направле- 





45 4-5 0 05 1 15 15 410,5 0 05 1 15 | 4,5 4-05 0 0,5 1 15 


Рис. 2. Решение внешней граничной задачи (4) для круга с г - 2 при ад - С ДР 
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ний градиента бі з 0,8315, бо - 0,5556 (как любое постоянное направление, зто не обеспе- 
чиваєт вьшолнения условия (12)) иллюстрируєт рис. З (М - 8 итераций). Точное решение 





цап У (хі ш 1) - фа М 1) - комплексно-сопряженная (й - и Ч їу - - голоморфная функ- 
ция) с ним функция у по соотношению Коши -- Римана также удовлетворяют зтой задаче. 
При б итерации (5) стремятся не к точному решению и, а к комплексно-сопряженной 
функции у. 

Вьщчисление аномалий. На практике отьскивают не потенциал у - и - у силь тяЯжести, а 
потенциал и аномалиий сильг тяжести (потенциал нормальной сильт тяжести у известен). В ли- 
нейньжх приближениях граничньх задач грави- и магнитометрии переход к аномалий не из- 
меняет не только уравнения, но и вид оператора граничньх условий |Интерпретация ..., 
1992). Необходимо лишь в качестве граничной функции м (5) взять аномалию силь тяжести. 

В нелинейньх граничньх задачах переход к аномалий изменяєт левую часть граничньх 
условий (З), (4), так как операция вьшчисления модуля градиента не дистрибутивна: | гад М | я 
т | ягад И | з | ягад У | . Если граничная функция м (5) - - модуль градиента потенциала у/ СИЛль 
тяжести, для потенциала аномалий и получаем нелинейную граничную задачу: 


3 
ди 
Ац-о, У, ах: ча, дхі дхі - у? (5)- | вайду| 7, побу» | (З 


і 21 і зі 


граничнье значения которой при нулевом потенциале нормальной сильг тяжести совпадают 
с граничньми условиями (3). Решают зту задачу одним из последовательньхх приближений 
(условия на го опущень(): 


ду (ди, | дЧачі 


2 з 
ої дХі| ді дХі | ага |? 


Мазі з 0, зу(5) 5 зно 








| агад (у ни )| у (5) 








диу/дх, 5-7 





/ 
РЕА 5 8 / 
0 04 08  -08 -04 0 04 08 





58 04 0 04 08 08 04 


Рис. 3. Решение внутренней граничной задачи (4) для единичного круга. 
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- ду | дип | Чан 
У 12 Б 
ді ді ді 














зі | вгаду |" 
Алі 70, УА (14) 

м (5) ч (5) 

5 

3 
- ду , ди дані 

ди ек, п БАЧ) -З В 4 2 
ші ря дХі ах дХі МРАУ у втаду 

Й 5 


Определение численной скорости сходимости итераций (14) имеет практический инте- 
рес, так как на каждой итерации решаєтся граничная задача и слабосходящиеся процес - 
сь существенно ухудшают зксплуатацию алгоритмов. Сходимость (14) вьічислена на модели 

1 








У - Ма , м (5) - | ягад (0 - У )| на границе 
2 2 2 
ху хо нХ о хо (хз я й)? 
5 (плоскость Хз - 5 на 36 узлах коллокации в точках, Х;| - - 25 -10К, Х, - - 25410), Х3 7-9, К, 


7-1,2,3, 4, 5) для 36 фундаментальньж решений, полюса которбьіх находятся под узлами 
коллокациий на плоскости Хз - 1) до вьшолнения условия тах | | ягад (У ЧИ) | -м (5) | «є. 
5 


При є- 10-93 первьшй итерационньй процесс завершился за три итерации, второй -- 36, по- 
следний -- 103 итерации (при є - 1031-- за 4, 78 и 146 соответственно). 

Для всех итераций (14) первой решена граничная задача с производной вдоль Хз на гра- 
нице. Вспомогательнье точки для фундаментальньтх решений взять на плоскости Х3 - 1. 
Приближение решения и его производньх из всех итераций хорошо совпадают между собой 
(при є - 103 с точностью 10 2 ,приє- 107" с точностью 1071; но далеки от точньх. Их сред- 
неквадратичньєе значения в узлах коллокации (1-й столбец) и в сАреа в между ними (2-й 
столбец) даєт табл. 3, а соответствующие погрешности для разньткх є- 10 Каре 1,2,..3)ми М -- 
табл. 4. 

Последняя таблица показьгваєт, что увеличение точности аппроксимации граничнькх усло- 
вий в узловьїх точках коллокации без увеличения количества узловьгх точек (увеличения чис- 
ла функций в разложений) не влияет на погрешность (начиная с К - 2). Зто говорит о нели- 
нейном поведениим погрешностей решения граничньх задач с производньтми и побуждаєт ис- 
следовать погрешности решений каждой конкретной задачи. 

Вьводьгю. Трансформацию значений поля силь тяжести во внешнем пространстве можно 
осуществить путем решения граничной задачи (3) для уравнения сильг тяжести (2), конкрет- 
ньй вид которого зависит от избранной модели Земли. В плоском случаєе зта операция сво- 
дится к решению граничной задачи (4) по итерациям (5). Точность ее решения существенно 
зависит от мерьг обусловленности задачи, приближений направления искомого градиента, 
параметров численного метода |Алексидзе, 1978| (числа и положения фундаментальньх ре- 


Таблица 3. Приближеннье решения нелинейной граничной задачи (14) 
на плоскости 
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Таблица 4. Погрешности решений нелинейной граничной задачи (14) 
на плоскости при разной точности 








шений), направления внешней нормали, вдоль которой вьтчисляют производнь6е в граничньх 
условиях (вьшолнение условия | 12)). Невзирая на зто, ее решение на простьгх моделях при 
определенньгх ограничениях успешно. Основнье источники неоднозначности решений зада- 
чи -- конечное вертикальное разрешение решений, неадекватньгй вьбор модели задачи и на- 
чальньжх приближений решения. 

Благодарю академика В. И. Старостенко за ценнье замечания касательно рукописи статьи. 
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